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The objective of this study is to verify the Symmetrical 
Number System (SNS) undersampling receiver architecture using 
software and investigate implementation issues using Digital 
Signal Processing (DSP) hardware. In the software design, a 
MATLAB program is written to determine a single sinusoidal 
input frequency using this receiver architecture. Each 
channel of the SNS undersampling receiver consists of a low 
speed ADC, a discrete Fourier transform followed by a 
constant threshold device to detect the signal's frequency 
bin. The detected frequency bins are then recombined in a 
SNS- to-decimal algorithm to recover the frequency of the 
signal . Error rate performance in a Gaussian noise 
environment at the input stage is evaluated. In the hardware 
design, a sinusoidal waveform is digitized, discrete Fourier 
transformed and converted from the SNS format to a decimal 
value using a single channel digital signal processor. 
Implementation difficulties and design issues are discussed. 
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I . INTRODUCTION 



A . UNDERS AMPL ING 

The digitization of a signal is usually governed by the 
Nyquist theorem where the sampling frequency is at least 
twice the signal bandwidth. The Nyquist theorem however, 
places a limitation only on the information that can be 
derived from a single set of digitized data [Ref. 1] . If the 
sampling frequency is less than twice the bandwidth of the 
signal being digitized, aliasing and consequently ambiguities 
occur. With additional information however, ambiguous 
frequency components due to undersampling may be resolved. 
Such information may come from, for example, trial sampling 
periods. Rader [Ref. 2] described how trial sampling periods 
can be used to recover periodic signals. The trial sampling 
period which yields the waveform of smallest variation is 
considered to be the correct period and the resulting 
waveform the correct waveform. 

Pace, Leino and Styer [Ref. 3] examined the relationship 
between the Discrete Fourier Transform (DFT) and the 
Symmetrical Number System (SNS) as a means of resolving 
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single frequency undersampling aliases. They, showed that the 
DFT encodes the frequency information of a signal in a format 
that is in the same form as the SNS . In addition, they 
proved analytically that aliases resulting from undersampling 
a single- frequency signal could be resolved using 2 or more 
channels. Each channel in a SNS undersampling receiver 
contains a low speed ADC, a DFT and a threshold device to 
detect the input signal bin number in the frequency domain. 
The bin numbers from each channel are then recombined to 
resolve the signal's frequency. 

B. PRINCIPAL CONTRIBUTIONS 

First, this thesis verifies the SNS undersampling theory 
advanced by Pace, Leino and Styer [Ref. 3] . An algorithm is 
written and coded in MATLAB to prove the methodology and to 
show that the frequency of an undersampled signal can be 
accurately measured. The algorithm is also simulated in a 
Gaussian noise environment. Error rates for the different 
noise levels are obtained as a function of the signal to 
noise ratio. Since the Fast Fourier Transform (FFT) is not 
suitable for computing DFTs in this application, alternative 
methods are suggested for real-time applications. 
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Second, possible hardware implementation problems are 
investigated based on a Digital Signal Processing (DSP) 
platform. Several problems were encountered: the need for 
stable sampling frequencies, large memories and alternative 
methods for computing DFT for fast response time. 
Integration into future EW receivers must take these factors 
into consideration. 

Undersampling offers several advantages [Ref. 4] . It 
allows the resolution of very high frequencies in EW 
receivers using low speed ADCs. This is especially so if 
several SNS channels are used. In particular, the use of 
undersampling in the design of receivers will reduce their 
cost and complexity. 

C. THESIS ORGANIZATION 

4 

In Chapter II, the relationship between the SNS and the 
digital frequency domain as mapped by the DFT is examined as 
a means of resolving single -frequency under sampling 
ambiguities. It shows how the frequency of a signal that is 
undersampled at two different sampling frequencies (two- 
channel) can be determined. In order to use lower sampling 
frequencies, the two- channel case can be extended to three or 



3 



more channels. 



In particular the three-channel case is 



discussed. 

In Chapter III, algorithms for the two- channel and 
three-channel receivers are developed and coded in MATLAB to 
measure the frequency of an incoming signal . Each section of 
the software is explained in detail. Results are obtained 
based on different Gaussian noise levels. 

A feasibility study/design for the two-channel case is 
carried out in Chapter IV using a DSP development kit . The 
suitability of using a DSP platform and its associated 
problems are discussed. 

Chapter V states some conclusions and recommendations 
for future research. 
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II. BACKGROUND INFORMATION 



A . INTRODUCTION 

Digitization of a signal is usually governed by the 
Nyquist criterion when the input signal is bandlimited to 0 < 
f < f s /2 where f s is the sampling frequency. For higher 
frequencies (i.e. f > f e /2) , the process of undersampling 
gives rise to ambiguities. However, with additional 
information (or channels) , the frequency components f > f s /2 
can be resolved. 

Pace, Ramamoorthy and Styer [Ref. 5] showed that the 
discrete Fourier transform (DFT) naturally encodes the 
frequency information of a signal in the same format as the 
symmetrical number system (SNS) . Consequently, aliases from 
undersampling can be resolved using this method. The theory 
set forth is elaborated in [Ref. 3] . 

B. DISCRETE FOURIER TRANSFORM (DFT) 

Since all signals consist of sinusoids, for simplicity, 
a single frequency sinusoidal waveform is used for analysis. 
Assume the sinusoidal signal is 

x(t) = 2 cos to t 
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( 1 ) 



and after sampling 

x(n) = 2cos© n. 



The DFT of x(n) is given by [Ref. 6] : 



( 2 ) 



X(k) = ^x(n)e" j(2mk/N) k = 0,1,... ,N - 1 . 

n=0 



(3) 



Applying the DFT to x(n) results in a discrete spectrum where 
|X(k)| 2 is the energy contained in the signal at each digital 
frequency co=27ik/N. The spectrum X(k) has N indices with the 
digital frequency of each index given by: 



1 ^(N/2) ^(N/2 + l) _ (N-2) (N-l) 

0,271— ,...,2k - yliz — — — — ,271 — — — 

N N N N N 



for N even 



(4) 



and 



■ ft , 1 o_(N-l)/2 (N + l)/2 

t; r: , 

N N N 






N 



N 



for N odd. 



(5) 
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The analog frequency corresponding to each index is obtained 
by multiplying each value by f B . Since signals with digital 

frequencies in the range n < © < 2tz are indistinguishable 

✓ 

from signals with digital frequencies 0 < © < k, the digital 
frequency of each index can also be written as: 



N N N N N 



for N even 



and 



( 6 ) 



' 1 „ In / 2 J |_N/2| . 2 . 1 

0,27c — ,...,271— -,2n- -,...^2k — ,27c — 

N N N N N 



for N odd. 



(7) 



where LxJ is the floor function and represents the greatest 
integer less than or equal to x. Thus the spectrum X(k) 

4 

resolves into N integer indices and incoming signals will map 
into unique bins : 



„ , N N , . , 



for N even, 



( 8 ) 





N 


N 




0 , 1 ,..., 


_2 j 


_T_ 


,..., 2,1 



for N odd. 



( 9 ) 
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For example, for N = 5 (f B = 5 Hz and the sampling duration 
T x is 1 second) , the output bins after the DFT are [0122 
1] for input frequencies of [01234] Hz. These DFT bins 
are repeated for higher frequencies as illustrated in Figure 
1. In this figure the abscissa corresponds to the incoming 
frequency and the ordinate corresponds to the bin into which 
the signal is resolved. 
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Figure 1: DFT bin mapping for input frequencies f= 

0 to 10 for N = 5 (f s = 5 Hz sampling for 1 second) . 
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C. THE SYMMETRICAL NUMBER SYSTEM (SNS) 

The SNS is composed of a number of pairwise relatively 
prime (PRP) moduli. The integers within each SNS modulus 
however, are derived from a symmetrically folded waveform. 
The symmetrically folded waveform corresponding to each SNS 
PRP moduli (mi), has a folding period equal to the modulus. 
The integer values within each SNS modulus are derived from a 
mid-level quantization of the symmetrical folding waveform. 
The formal definition of a symmetrical residue is given 
below: 

Definition: For an integer h such that 0 < h < m 

x h = min {h, m - h} 

( 10 ) 

If this function is extended periodically with period m, 
that is , 

X h + nm = X h 

( 11 ) 

where n € {0, ±1, ±2 , .„} then x h is called a symmetrical 

residue of (h+nm) modulo m. For m even, let x be the row 
vector 
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X = 



^ , mm „ ^ „ 

~~ •••> 2, 1 



( 12 ) 



For m odd, let x be the row vector 







m 


m 




x = 


0,1,..., 


y] 


_T_ 


,.-,2,1 



(13) 



where LxJ again represents the floor function resulting 
in the greatest integer less than or equal to x. These 
two vectors consist of the symmetrical remainder 
elements x h , 0 < h < m. 



D. RELATIONSHIP BETWEEN DFT AND SNS 

From the above, it is obvious that the DFT maps real 
signals naturally into the SNS. That is, in Section C, if we 
let the modulus m represent the sampling frequency multiplied 
by the sampling time (i.e., f a T i ) , then equations (12) and 
(13) are in the same form as equations (8) and (9) where N= 
f e T x . Thus the SNS provides a convenient framework for 
undersamping signal analysis. 

Table 1 displays the input frequencies and the resulting 
DFT bins for sampling frequencies 5 Hz and 6 Hz respectively. 
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Input Frequency 


DFT Bins 


f 


f s = 5 Hz | f s = 


6 Hz 


o 


0 ! 


o 


1 


i ! 

t ... 


1 


! 2 


2 ! 


2 


3 


2 


3 


4 


1 


2 


5 


0 ! 

1 


1 j 


P 6 

: 


1 


0 


7 


2 i 


1 


8 


2 


2 



Table 1: Input Frequency and Resulting DFT Bins for 2 

Channel Example . 

The frequencies are resolved as described in equations (12) 
and (13) . By considering two or more channels, it is 
possible to unambiguously resolve the signal frequencies in 
the dynamic range determined by the SNS. One method is to 
devise a look-up table similar to that shown in Table 1 . 
However this method is inefficient for high frequencies; 
large memories are required. An alternative method is 
described below: 

Suppose there are r channels and the incoming frequency 
is within the dynamic range of the system. To carry out the 
SNS-to-decimal conversion, we need to solve f = a A (mod 11 ^) 
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for i 



1,2,... r, where a ± is the corresponding detected DFT 



bin for each n^. The Chinese Remainder Theorem states that 
there is a unique solution modulo M=m 1 *m 2 . ..* 1 ^. A standard 
method of solution is to find integers b ± such that M*b i /m i = 
1 (mod m A ) where i = l,2,....r in which case the solution is f 
= M*b 1 *a 1 /m 1 + M*b 2 *a 2 /m 2 + . . . + M*b r *a r /m r (mod M) . In Sections 
F and G below, examples are given to illustrate this 
calculation. 

E. DYNAMIC RANGE OF THE SNS 

Let m 1( . . . , m,. be r pairwise relatively prime moduli, 

then the dynamic range, D (0:D-1) of a SNS system is given as 
follows : 

• If all the moduli are odd, then the dynamic range of 



the system is 




(14) 



where j ranges from 1 to r-1 and m i2 ,m l3 ...m^ range over 



all permutations of {l, 2 , 3 , ..., r} . For example, for a 



two-channel case with m x = 5, = 7, 
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„ . I’m. m, 

D = mini — M — - 
1 2 2 . 



or D = 6. 



For a three-channel case with m^, = 3, to, = 5, m 3 = 7, 



D = ^-min{m, + m 2 m 3 ,m 2 + 111^3,1113 + m,m 2 } 



or D = 22. 



If one of the moduli (m x ) is even, then the dynamic 



range of the system is 



D = mini— F 



l 2 >= 



n^+rK 



i=+i 



(15) 

where j ranges from 1 to r-1 and m,,m i3 ...m lr range over 
all permutations of {2,3,...,r}. For example, for a 
two-channel case with m x = 6 to, = 5, 



K- 



D = min< — L +m 

l 2 



or D = 8. 
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For a three -channel case with tr^ = 8, = 5, nij = 7, 



f m, m 1 m, 

D = rninj + m ; m 3 , -^-m 2 + m 3 , -^-m 3 + m 

or D= 27. 

Clearly, the dynamic range of an SNS system with one 
even modulus is superior to that using all odd moduli. 
Moreover, the greater the number of channels, the greater the 
dynamic range . 

F. THE TWO-CHANNEL CASE 

Figure 2 shows the block diagram of a two- channel 
receiver architecture to determine a single frequency f. In 
this architecture the ADC sampling frequencies f sl and f s2 are 
relatively prime and T 3 = 1. The DFT outputs are thresholded 

t 

to detect the frequency bins of the signal. The detected 
frequency bins a x and a 2 are then used by the SNS-to-decimal 
algorithm to determine the frequency of the input signal . 
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frequency bins a x and a 2 are then used by the SNS- to -decimal 
algorithm to determine the frequency of the input signal. 



ADC 

fsl 



N|/ 



Window 
Funct ion 



DFT 



Bin 

Detector 



!\nti -Aliasing 
Filter 



ADC 

fs2 



Window 

Function 



DFT 



Bin 

Detector 



a2 




Figure 2 : Block Diagram of a Two Channel Receiver 
Architecture . 

Let m x = f sl and = f B2 and suppose that the incoming 
frequency is within the dynamic range of the system. From 
Section D, we need to solve f = a^mod m 2 ) and f = a 2 (mod n^) . 
The two congruence equations, f = a^mod m 2 ) and f = a 2 (mod 
trij) are solvable only if the greatest common divisor of m x 
and divides (a 2 - a 2 ) , a generalization of the Chinese 

Remainder Theorem [Ref. 7] . To solve for f, the diophantine 
equation 

p^+q+nu, = (a 2 - aj (16) 
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must be solved for p and f is then calculated from the 
equation 

✓ 

f = a 1 +p*m 1 . (17) 

The code for this algorithm is shown in Appendix A. 

For example, for sampling frequencies 5 and 6, m x and rr^ 
have values of 5 and 6 respectively (T x = 1) . If the signal 
is resolved into bins a 2 (=2) and a 2 (=1) after the DFT, p 
is found to have a value of 1 and q is found have a value of 
-1. Thus, the input frequency from (17) is 2+1*5 = 7. This 
can also be verified as shown in Table 1. 

G. THE THREE-CHANNEL CASE 

Figure 3 shows the block diagram of a three -channel 
receiver architecture to determine a single frequency f . 

4 

Similar to the two-channel case, the ADC sampling frequencies 
f.i, f B2 , and f s3 are pairwise relatively prime and T x = 1. The 
DFT outputs are thresholded to detect the frequency bins of 
the signal. The frequency bins a 2 , a 2 and a 3 are then used by 
the SNS-to-decimal algorithm to determine the frequency of 
the input signal. 
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Figure 3 : Block Diagram of a Three Channel 
Receiver Architecture. 

In the three -channel solution, let = f sl , = f s2 and 
trig = f e3 and suppose that the incoming frequency is within the 
dynamic range of the system. We need to solve f = a 1 (mod m 2 ) 
and f = a 2 (mod mj and f s a 3 (mod to,). Using the Chinese 
Remainder Theorem and the Euclidean algorithm, the method of 
solution is to find integers b ± such that M*b i /m i s 1 (mod raj 
where i = 1,2, and 3 and M = m 2 * * rr^. The solution is 

then f = ± M*b 1 *a 1 /m 1 ± M*b 2 *a 2 /m 2 + M*b 3 *a 3 /m 3 (mod M) where f 
is the frequency which falls within the dynamic range D of 
the system. 
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For example, let m x = 5, = 6 and rr^ = 7, so that M = 

210 and D = 22 . Suppose that the signal is resolved into 

bins a 3 (= 1) , a 2 (= 2) and a 3 (=2) after the DFT . For the 

✓ 

three-channel case the b^ values must be found. Here, b x b 2 
and b 3 are found to be -2, -1, and -3 respectively. Thus f = 
± 210 (-2) (1) /5 + 210 ( -1) (2) /6 ± 210(-3)(2)/7 mod(210) and we 
must choose the solution that falls within the SNS dynamic 
range D = 22 [0:21] . The correct combination f = 84 -70 + 

180 mod (210) = 194 mod (210) . Although 194 is out of the 

dynamic range, 210 - 194 = 16 is in the dynamic range so that 
f = 16 is the correct frequency. 



H. NOISE CONSIDERATIONS 

For a sinusoidal waveform, the Signal to Noise Ratio 
(SNR) is defined as 



P 

SNR = — =- 
2a 2 . 



(18) 



where P is the power of the signal and a 2 is the noise power. 
Assuming a signal power of one, the noise power and amplitude 
are given by 
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1 



2 SNR 



( 19 ) 



1 

72 SNR. 



( 20 ) 



This a is multiplied by a normally distributed random number 
sequence of zero mean and unit variance and added to the 
input signal as noise. The simulation results are given in 
Chapter III. 
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III. SOFTWARE DESIGN AND RESULTS 



A. TWO-CHANNEL ALGORITHM 

The two-channel case was described in Chapter II. An 
algorithm was constructed based on Figure 2 . The software 
given in Appendix A can be divided into the following 
sections : 

• Initialization. This section obtains all the 

parameters (number of iterations, input frequency, 
sampling frequencies, quantization levels) required. 

• Iteration loop. This section consists of a loop 
(with an initial count of zero) to count the number 
of errors . 

• Creation of Waveform. Based on the input frequency, 
a sinusoidal waveform is created with noise added. ' 

• Sampling and Quantization. The waveform is then 
sampled at two different frequencies and quantized 
using a 14 -bit ADC. 

• Windowing. A rectangular window operation of width N 
= f B * Ti = f B (the total sampling/ integration time is 
taken to be one) is carried out. 
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• DFT Operation. A DFT is then carried out on each 
sample, taking only the first half of the DFT output. 
The formula used for the DFT process is a simple pair 
of nested loops . 

• Bin Detection. A non-adaptive (constant) threshold 
bin detector is then used to find the bin with the 
maximum value for each DFT output . 

• SNS-to-Decimal Algorithm. The SNS- to -decimal 

algorithm as described in Chapter II is then used to 
calculate the incoming frequency. 

A flow diagram of this algorithm is illustrated in 

Figure 4. The MATLAB code can be found in Appendix A. 
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Figure 4 : Two Channel Algorithm 
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B. TESTING OF TWO-CHANNEL SYSTEM 

To test the two-channel case (sinusoidal signal without 
noise) , the program is run with the following input and 
sampling frequencies shown in Table 2. 



f 




f. 2 


Dynamic Range 


Remarks 


7 


5 


8 


0:8 


Low input frequency 


100 


97 


98 


0:145 


Consecutive sampling 
frequencies 


1040 


547 


1200 


0:1146 


Sampling frequencies 
far apart 


12125 


12671 


12919 


0:12794 


High input frequency 



Table 2: Tested Input and Sampling Frequencies. 

For example, with input signal frequency at 7 Hz as 
shown in Figure 5, the sampled signals at 5 Hz and at 8 Hz 
are shown in Figures 6 and 7 respectively. The DFT output 
for the two samples are shown in Figures 8 and 9 . The 
resultant bins of the first halves of Figure 8 and 9 are then 
supplied to the SNS-to-decimal algorithm to be converted to 
the input frequency of 7 Hz. 
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Amplitude 




Figure 5: Input signal with frequency of 7 Hz. 
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Figure 6 : Sampled signal at frequency 5 Hz . 
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Amplitude 




Figure 7: Sampled signal at frequency 8 Hz. 
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Figure 8 : DFT output with f sl=5 Hz . 
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^ 6 - 




magnitude 




Figure 9: DFT output with fs2=8 Hz. 
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It is found that if one of the sampling, frequencies was 
the same as the input frequency, the algorithm failed. This 
is because the resulting samples due to the same sampling 
frequency will consist of zeros. This problem can be solved 
by using at least two sets of sampling frequencies. Apart 
from this, the algorithm works well in this noise-free (high 
signal -to -noise ratio) environment. 

C. SIMULATION PARAMETERS FOR TWO-CHANNEL CASE 

To obtain the error rates in a noisy environment, the 
two-channel software is run with the following parameters: 

• Number of iterations, num = 10000 

• Signal to Noise Ratio, SNRDB = -30 to 30 dB 

• ADC resolution, bit = 14 

• Input and sampling frequencies as shown in Table 3 . 



f 




f s2 


9 


10 


11 


90 


91 


92 


900 


901 


902 


9000 


9001 


9002 



Table 3: Input and Sampling Frequencies. 

D. RESULTS FOR TWO -CHANNEL CASE 

The results obtained are shown in Figure 10. 
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Figure 10 : Error Rates vs . SNR for two-channel system 
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The following observations are made: 



• As expected, the error rates improve as the SNR 
increases. A tradeoff between SNR and error rate is 
required. 

• Improvements in error rates were obtained when higher 

frequencies were used. This is because at higher 

frequencies, higher sampling frequencies are 
required. This leads to a higher N-point DFT (higher 
gain) which is less affected by noise. 

• However, at higher frequencies, the time taken to 

compute the DFT was much longer. To reduce the time 

taken, the following methods can be implemented: 

• If N is highly composite (factorable into powers of 
many small prime factors, preferably primes < 10) , 
use a "mixed-radix" FFT implementation. 

• If N is prime, or contains very large prime 

factors, use the "chirp-z" transform. 

• Use three or more channels in the receiver. A 

1 

three-channel receiver has a higher dynamic range 
for the same magnitude of sampling frequencies . 
For example, a two-channel receiver with sampling 



32 



frequencies 6 and 7 has a dynamic range of [0:9] 
while a three -channel receiver with sampling 
frequencies of 5, 6 and 7 has a dynamic range of 

[ 0 : 21 ] . 

THREE -CHANNEL ALGORITHM 

The three -channel algorithm is similar to the two- 
channel algorithm as shown in Figure 11. The MATLAB code 
can be found in Appendix A. 




Figure 11 : Three Channel Algorithm 
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F. TESTING OF THREE CHANNEL ALGORITHM 

To test the three-channel case, the program is run with 
some of the following input and sampling frequencies in Table 
4. 



f 


fsx 






Dynamic 

Range 


Remarks 


13 


5 


6 


7 


0:21 


Low input frequency 


100 


17 


18 


19 


0:171 


Consecutive 

sampling 

frequencies 


1040 


17 


91 


919 


0:1232 


Sampling 
frequencies far 
apart 


12125 


90 


929 


937 


0:42741 


High input 
frequency 



Table 4: Tested Input and Sampling Frequencies. 



Apart from the anomaly discussed in the two- channel 

* 

case, the algorithm works well in this noise- free (high 
signal -to-noise ratio) environment. 

G. SIMULATION PARAMETERS FOR THREE -CHANNEL CASE 

To obtain the error rates in a noisy environment, the 
three-channel software is run with the following parameters-. 

• Number of iterations, num = 10000 
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Signal to Noise Ratio, SNRDB = -30 to 30 dB 



• ADC resolution, bit = 14 

• Input and sampling frequencies as shown in Table 5 . 



f 


fsl 


f s2 


f 8 3 


9 


5 


7 


11 


90 


13 


14 


17 


900 


41 


42 


43 


9000 


141 


142 


143 



Table 5: Input and Sampling frequencies. 

H. RESULTS FOR THREE - CHANNEL CASE 

The results obtained are shown in Figure 12 . 
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Figure 12 : Error Rates vs . 



SNR for three-channel system 
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Comparing the two- channel and three-channel cases, the 
following observations can be made: 

• The three-channel system is much faster than the two- 
channel system since the DFTs required are smaller 
due to the smaller sampling frequencies. 

• However the results for the two-channel system with 
noise are better. For example to achieve a 
relatively error-free system for a frequency of 9000 
Hz, the two-channel case requires only -22 dB. 
However, the three-channel case requires at least -4 
dB. 
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IV. HARDWARE DESIGN AND FINDINGS 



A . INTRODUCTION 

In the last chapter, the instantaneous measurement of 
frequency using the SNS- to -decimal algorithm was verified. 
There is a need to investigate the implementation of the 
algorithm in hardware. Digital Signal Processing (DSP) 
hardware was selected for the following reasons: 

• A major part of the algorithm is the processing of 
DFTs which is a digital signal processing task well 
suited to be carried out by DSP hardware. 

• DSP hardware provides a fast way to implement the 

algorithm. The DSP development kit is easy to learn, 
program and simulate. It is ideal for this 

application to investigate hardware problems and 
limitations . 

• Cost consideration: the development kit plus tools 

cost $1500; 

• EW receivers are likely to incorporate DSP hardware. 
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B. TI TMS320C54X DSP DEVELOPMENT KIT 

The TMS32054C54X DSKplus [Ref. 8-12] is a low cost DSP 
starter kit that gives a designer a working knowledge of DSP 
code to build DSP based systems . The development kit 
contains a stand-alone application board that can be 
connected to the PC. It executes code in real time at 40 
MIPS while the Windows-based debugger analyzes it line-by- 
line, displaying internal DSP register information in 
multiple windows and in real time. It has an Analog 
Interface Circuit for the input of signals. The board's 
communication interface enables the creation of C54x DSP code 
and host PC code. Moreover, the hardware enables the use of 
expansion slots for adding memory, peripherals such as 
interface logic, other DSPs etc. The developed code can 
eventually be loaded into a resident DSP processor, which may 
be part of a EW receiver architecture . Figure 13 shows a 
block diagram of the development kit . A more detailed 
description of the kit can be found in Appendix B. 
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Parallel Port to PC To Other Interface Ports 
Figure 13 : Block Diagram of DSP Hardware . 



C . SOFTWARE 

The software for the two-channel case described in 
Chapter II (Figure 2) is written using the DSP development 
kit. The software (found in Appendix C) is coded in 'C' 
language/assembly language and converted to the C54x assembly 
language (if required) prior to execution: 

• Firstappl . c/Firstapp2 . c . These two programs poll the 
input channel and sample the input signal at the two 
sampling frequencies respectively. 

• Hostappl.cpp/Hostapp2.cpp. These two programs 

display the samples of the signals based on the two 
sampling frequencies and save the data in text files. 
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• Main . c . 



This program reads the data, executes the 



DFT, obtains the largest values for the two channels 
and then carries out a SNS- to -decimal conversion. 

These programs were run individually and consecutively. 

D. TESTING AND RESULTS 

Using data generated by MATLAB, the main program was 
tested successfully in the development kit. The programs 

were then run with an input frequency of 126 Hz and sampling 
frequencies, 125 Hz and 128 Hz. Results obtained were 
intermittent i.e., correct results were not always obtained. 
A frequency counter and an oscilloscope were set up and it 
was found that the sampling frequencies were not stable. 
Testing with different frequencies did not improve the 
results . 

E . PROBLEMS 

Several problems were encountered during the 

investigation: 

• Stability of Sampling Frequencies. The development 
kit carries out frequency division of the master 
oscillator to obtain the sampling frequencies . 
Unfortunately, the crystal oscillator has a 
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resolution of 5-10 Hz. 



This is unacceptable as a 



shift of 1 Hz in the sampling frequency will cause 
erroneous results. Moreover, the fact that the 
sampling frequencies are factors of the oscillator 
frequency and that they need to be pairwise 
relatively prime severely limits the choice of 
frequencies . A possible solution is to obtain the 
sampling frequencies directly from stable signal 
sources . 

• DFT. For higher frequencies, the execution of the 
DFT takes a long time. Several solutions were 
suggested and discussed in the previous chapter. 

• Memories. Insufficient memory error messages were 
encountered when high frequencies were used. The 
same messages occurred when attempts were made to run 
the routines together. More memories and/or more 
efficient DFT algorithms are required. 
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V. CONCLUDING REMARKS 



The main contribution of this thesis is the verification 
of the relationship of the DFT to the SNS to resolve 
undersampling ambiguities and the investigation of hardware 
implementation issues using a DSP platform. Error rates for 
different SNR are also obtained. 

The use of undersampling technique using the SNS to 
measure frequency is a viable method to implement in a EW 
receiver architecture. However, the need for faster DFT 
computation and stable sampling frequencies must be taken 
into account before they can be considered for incorporation 
into EW receivers . There is also a trade-off between the 
number of channels and SNR. For faster response, a multi- 
channel case is recommended; but a higher SNR is required. - 
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APPENDIX A 



MATLAB CODE FOR SOFTWARE ALGORITHM 



% Thesis Project 
% 

%■ Two Channel Receiver 
% 

% Note: The sampling frequencies should be relatively prime 



clear all? 

% Initialization 



num=input ( ' Enter Number of iterations : ' ) ? % 
f=input ( 1 Enter Input Frequency: ' ) ; % 
fsl=input (' Enter Sampling Frequency 1:'); % 
f s2=input ( 'Enter Sampling Frequency 2: ' ) ? % 
fpl=fopen ( 1 c:\matlab\bin\thesis\result.dat ' # 1 at ' ) ? 



Number of iterations 
Frequency of signal 
Sampling frequency 1 
Sampling frequency 2 
% Store results 



% Quantization levels 



bit=14 ? 

qnlevel=2 A bit-l ; 
q=2/qnlevel; 



% No. of quantization levels 
% quantization size 



for SNRDB= -30:2:30 



% Set Signal to Noise Ratio 
% from -30 dB to 30 dB 



count =0 ; 



% Error Count 



for 



i=l :num 



SNR=10 (SNRDB/10) ; 
sigmasq=l/2/SNR; 

t= ( 0 : .001:1) ; 
sig=sin (2*pi*f *t ) ; 
tl=l/f Sl : l/f Sl : 1 ; 
noisel=sqrt (sigmasq) *randn(l, length 
ADCsigl=sin (2*pi*f *tl) +noisel; 
ADCsigl=f ix (ADCsigl/q) *q? 
t2=l/fs2:l/fs2:l; 

noise2=sqrt (sigmasq) *randn (1, length 
ADCsig2=sin (2*pi*f *t2) +noise2; 
ADCsig2=f ix (ADCsig2/q) *q? 



% Convert to non-dB units 
% Noise normalization assuming 
% signal power of 1 

% signal 
% first ADC 
tl) ) ; % noise 
% digitized signal 
% quantized signal 
% second ADC 
t2) ) ; % noise 
% digitized signal 
% quantized signal 



%f igure (1) 

%subplot (3 , 1, 1) , plot (t,sig) 
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%title (' Figure 1. Plot of signal') „ 

%xlabel ( 'Time' ) 

%ylabel ( 'Amplitude' ) 

% subplot (3,1,2) ,plot (tl , ADCsigl) 

%title ( ’ Figure 2. Plot of sampled signal (sampling frequency 1) 
plus noise 1 ) 

%xlabel ( ’ time’ ) 

Irylabel ( * magnitude 1 ) 

%subplot (3,1,3) ,plot (t2,ADCsig2) 

%title (’ Figure 3. Plot of sampled signal (sampling frequency 2) 
plus noise’) 

%xlabel ( * time ’ ) 

%ylabel ( * magnitude ’ ) 

% Window operation 
% Assume rectangular window 

winsizel=fsl; 
winsize2=f s2; 

winsigl=ADCsigl (1 twinsizel) ; 
winsig2=ADCsig2 (1 :winsize2) ; 

% DFT Operation 

DFTsigl=abs (f ft (winsigl , winsizel) ) ; 

DFTsig2=abs (ff t (winsig2 , winsize2) ) ; 

DFTsigla=DFTsigl (1 : length (DFTsigl) /2 +1); % Taking half of image 
DFTsig2a=DFTsig2 (1: length (DFTsig2) /2 +1); % Taking half of image 

%f igure (2) 

% Plot to locate position of maximum value 

% Note that due to MATLAB (which cannot have a zero index, the 
actual location is one less 
%subplot (2,1,1) , stem(DFTsigla) 

%title (’ Figure 1. DFT plot of signal with sampling frequency 1’) 
%xlabel ( ’ frequency bins ’ ) 

%ylabel ( ’ magnitude ’ ) 

%subplot (2,1,2) , stem(DFTsig2a) 

%title (’ Figure 2. DFT plot of signal with sampling frequency 2’) 
%xlabel ( ’ frequency bins ’ ) 
lylabel ( ’ magnitude ’ ) 

% bin detector 

[i ,yl] =max (DFTsigla) ; % yl, y2 are locations of max values 

[ j ,y2] =max (DFTsig2a) ; % Note that due to MATLAB, the 

% actual location is one less. 

al=yl-l ; 
a2=y2-l ; 

% SNS to Decimal Algorithm 



% size of window is fsl 
% size of window is fs2 
% windowed sampled signal 1 
% windowed sampled signal 2 
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% To solve for f ==ai (mod mi) (where "= = " indicates congruence and 

V mi are pairwise relatively prime) , the Chinese Remainder Theorem 
% states that there is a unique solution modulo M=ml*m2 . . . .mr . 

% A standard method of solution is to find integers bi such that 
% M*bi/mi==l (mod mi) where i=l,2,....r in which the solution is 

V f==M*bl*al/ml + M*b2*a2/m2 + . . . + M*br*ar/mr (mod M) 

% For a 2 channel case, i.e. i=l,2, 

% m2*bl == 1 (mod ml) 

% ml*b2 == 1 (mod m2) 

% f == al (mod ml) 

% f == a2 (mod m2) 

% f == m2*bl*al + ml*b2*a2 (mod ml*m2) 

% Given ml (sampling frequency 1) and m2 (sampling frequency 2), to 
% find bl and b2 , the congruence equation is transformed to a 
% diphantine equation and solved using the Euclidean algorithm: 

% m2*bl - ml*yl = 1 

% ml*b2 - m2*y2 = 1 

% The above two equations can be combined into 
% m2*bl - ml*b2 = 1 

% bl and b2 are solved by the function n lde.m" which is called by 
"glde .m M . 

% 

% f == al (mod ml) and f == a2 (mod m2) is solvable only if the 
% greatest common divisor of ml and m2 divides (a2 - al) . 

% To solve for f, r from the diophantine equation 
% r*ml+s*m2 = a2 - al must be solved. 

% r is obtained from "glde.m" and f is calculated by the 
%• equation f = al+r*ml 

idif f=a2-al ; 

r = glde (fsl , fs2, idiff ) ; 
freq=abs (al+r*fsl) ; 

% Count the number of correct results. 



if freq==f 

count =count+l ; 

end 

end 

error = 1 -count /num; 

% Write results to file 

xl=fprintf (fpl , ' %d %d I’d %d %d I’dXn' , f, fsl, fs2, SNRDB, num, error); 
plot (SNRDB, error, *y+') 

title ('Error Rate vs. Signal to Noise Ratio') 
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xlabel ( 1 SNR(dB) ' ) 
ylabel ( ' Error Rate % ’ ) 
hold on 
end 

f close (fpl) ; 

% To calculate the dynamic range 

if rem(fsl , 2) ==0 % To check whether fsl is even 

DR=f Sl/2 + fs2; 

elseif rem(fs2 / 2) ==0 % To check whether fs2 is even 

DR=f s2/2 + fsl; 

else 

DR=. 5* (fsl+fs2) ; % fsl and fs2 are odd numbers 

end 

o, 

o 

% Thesis Project 
% 

% Three Channel Receiver 
% 



clear all; 
close 

% Initialization 

num=input ( 1 Enter Number of iterations : 1 ) 
f= input ( ' Enter Input Frequency: f ) ; 
fsl=input (' Enter Sampling Frequency 1:*) 
fs2=input (' Enter Sampling Frequency 2:*) 
fs3=input (' Enter Sampling Frequency 3: f ) 



% Number of iterations 
% Frequency of signal 
% Sampling frequency 1 
% Sampling frequency 2 
% Sampling frequency 3 



fpl=fopen( ! c:\matlab\bin\thesis\result .dat 1 , 'at ' ) ; 

% Store results in file for later processing if required 

% Quantization levels 

% bit=input (' Enter ADC resolution:'); 
bit=14 ; 



qnlevel=2 *bit - 1 ; 
q=2/qnlevel ; 

for SNRDB= -30:2:30 



% No. of quantization levels 
% quantization size 

% Set Signal to Noise Ratio from -30 
% dB to 30 dB 



count=0; 
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for i=l:num 

SNR=10^ (SNRDB/10) ; %• Convert to non-dB units 

sigmasq=l/2/SNR; % Noise normalization assuming 

% signal power of 1 



tl=l/f si : l/fsl : 1; % 

noisel=sqrt (sigmasq) *randn (1 , length (tl) ) ; % 

ADCsigl=1000* (sin (2*pi*f *tl) +noisel) ; % 

ADCsigl=f ix (ADCsigl/q) *q; % 

t2=l/fs2: l/f s2 : 1 ; % 

noise2 = sqrt (sigmasq) *randn (1, length (t2) ) ; % 

ADCsig2=1000* (sin (2*pi*f *t2) +noise2) ; % 

ADCsig2=f ix (ADCsig2/q) *q; % 

t3=l/fs3:l/fs3 : 1; % 

noise3=sqrt (sigmasq) *randn (1, length (t3) ) ; % 

ADCsig3=1000* (sin (2*pi*f*t3) +noise3) ; % 

ADCsig3=f ix (ADCsig3/q) *q; 

%f igure (1) 

%subplot (3,1,1) ,plot (tl , ADCsigl (l;fsl) ) 

%title (' Figure 1. Plot of sampled signal (sampling frequency 
1) plus noise') 

%xlabel ( ' time ' ) 

%ylabel ( ' magnitude ' ) 

%subplot (3,1,2) , plot ( t2 , ADCsig2 ( 1 : f s2) ) 

%title (' Figure 2. Plot of sampled signal (sampling frequency 2) 
plus noise') 

%xlabel ( ' time ' ) 

%ylabel ( 1 magnitude ' ) 

%subplot (3,1,3) ,plot (t3 ,ADCsig3 (l:fs3) ) 

%title (' Figure 3. Plot of sampled signal (sampling frequency 3) 
plus noise') 

%xlabel ( ' time ' ) 

Irylabel ( 'magnitude' ) 

% Window operation 
% Assume rectangular window 



first ADC 
noise 

digitized signal 
quantized signal 

second ADC 
noise 

digitized signal 
quantized signal 

third ADC 
noise 

digitized signal 



winsizel=f si; 
winsize2=f s2 ; 
winsize3=f s3 ; 

winsigl=ADCsigl (1 rwinsizel) ; 
winsig2=ADCsig2 (1 :winsize2) ; 
winsig3=ADCsig3 (l:winsize3) ; 



% size of window is fsl 
% size of window is fs2 
% size of window is fs3 
% windowed sampled signal 1 
% windowed sampled signal 2 
% windowed sampled signal 3 



% DFT Operation 

DFTsigl=abs (f f t (winsigl, winsizel) ) ; 
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DFTsig2=abs (ff t (winsig2 , winsize2) ) ; ^ 

DFTsig3=abs (fft (winsig3 , winsize3) ) ; 

DFTsigla=DFTsigl (1 : length (DFTsigl) /2 +1); % Taking half the image 

DFTsig2a=DFTsig2 ( 1 : length (DFTsig2) /2 +1); % Taking half the image 

DFTsig3a=DFTsig3 ( 1 : length (DFTsig3 ) /2 +1); % Taking half the image 

%f igure (2) 

%Plot to locate position of maximum value 

%Note that due to MATLAB (which cannot have a zero index, the 
%actual location is one less 

%subplot (3,1,1), stem(DFTsigla) 

%title (' Figure 1. DFT plot of signal with sampling frequency l 1 ) 
%xlabel ( 1 frequency bins 1 ) 

%y label ( 1 magnitude 1 ) 

%subplot (3,1,2) , stem(DFTsig2a) 

%title (' Figure 2. DFT plot of signal with sampling frequency 2') 
%xlabel ( ' frequency bins ’ ) 
lylabel ( 1 magnitude 1 ) 

%subplot (3,1,3) , stem(DFTsig3a) 

%title (' Figure 3. DFT plot of signal with sampling frequency 3') 
%xlabel ( ' frequency bins ’ ) 

%ylabel ( ' magnitude 1 ) 

% bin detector 

[i ,yl] =max (DFTsigla) ; % yl, y2 and y3 are the locations of 

% maximum values 

[ j ,y2] =max (DFTsig2a) ; % Note that due to MATLAB, the actual 

% location is one less. 

[k,y3] =max (DFTsig3a) ; 



al=yl- 1 ; 
a2=y2-l ; 
a3=y3 -1 ; 



% SNS to Decimal Algorithm 

bl=lde (fs2*fs3,fsl) ; 
b2=lde (f sl*f s3 , f s2 ) ; 
b3=lde (fsl*f s2 , f s3) ; 

cl=bl*fs2*fs3 ; 
c2=b2*f sl*f S3 ; 
c3=b3*f sl*f s2 ; 

f reqmat= [al*cl+a2*c2+a3*c3 ;al*cl+a2*c2-a3*c3 ; al*cl-a2*c2+a3*c3 ; 

al*cl-a2*c2-a3*c3 ; -al*cl+a2*c2+a3*c3 ; -al*cl+a2*c2-a3*c3 ; 
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-al*cl-a2*c2+a3*c3 ; -al*cl-a2*c2-a3*c3] ; 



freqmat=rem (freqmat, fsl*fs2*fs3) ; 

for i=l : 8 

if (freqmat (i) <0) 

fregmat (i) =freqmat (i) +f sl*f s2*fs3 ; 

end 

end 

freq=min(abs (freqmat) ) ; 

% Count the number of correct results. 

if f req-=f 

cou n t = cou nt + 1 ; 

end 

end 

error = 1-count/num; 

% Write results to file 

xl=fprintf (fpl, ? %d %d %d %d %d %d %d\n' , f, fsl, fs2, fs3, SNRDB, num, 
error) ; 

plot (SNRDB, error, *y+ ? ) 

hold on 

end 



fclose (fpl) ; 



% To calculate the dynamic range 



if rem(fsl,2) ==0 

x=[fsl/2 + fs2*fs3; fsl*fs2/2 
elseif rem(fs2, 2) ==0 

x= [f s2/2 + fsl*fs3; fsl*fs2/2 
elseif rem(fs3, 2) «0 

x= [f s3/2 + f s2* f si ; f s3*f S2/2 

else 



% To check whether fsl is even 
+ fs3; fsl*fs3 + fs2] ; 

% To check whether fs2 is even 
+ fs3; fs2*fs3 + fsl] ; 

% To check whether fs3 is even 
+ fsl; fsl*fs3 + fs2] ; 

Sr fsl,fs2 and fs3 are odd 



x=l/2* [fsl + fs2*fs3; fs2 + fsl*fs3; fs3 + fs2*fsl] ; 



end 

DR=min(x) ; 
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This function solves the general linear diophantine equation 
m2*bl - ml*b2 = k and returns the value bl 



function a=glde (ml , m2 , k) 

%• Calls function "lde" to calculate bl, b2 and na 



[bl,b2,na] =lde (ml, m2) ; 



% To check whether the equation is solvable. 

% na must be a factor of k for the equation to be solvable. 

mult=k/na; 



if (k-mult*na) ==0 % Equation is solvable 

bl=bl*mult; % These new values solve the diophantine equation 
b2=b2*mult ; 



mtest=bl; % To check whether bl and b2 are the least values 

mdl=ml/na; % that satisfies the diophantine equation 

md2=m2/na; 

mx=bl ; 

mx=mx+md2 ; 



while (abs (mx) -abs (bl) ) <0 
bl=mx; 
b2 =b2 -mdl ; 
mx=mx+md2 ; 

end 



if (mtest-bl ) ==0 
mx=bl; 
mx=mx-md2; 

while (abs (mx) -abs (bl) ) <0 
bl=mx ; 
b2=b2 +mdl; 
mx=mx-md2 ; 

end 

end 

end 



a=bl ; 



56 



% This function solves the linear diophantine equation 
% ml*bl + m2*b2 = na where ml and m2 are the sampling frequencies 
% and na is the greatest common divisor 
% and returns the value bl, b2 and na 
% 

% ml and m2 are assumed positive 

/ 

function [bl , b2 ,na] =lde (ml, m2) 

% Initialize bol, bo2, bl and b2 



bol=l ; 
bo2=0 ; 
bl = 0 ; 
b2 = l ; 

% Place ml and m2 in ma (dividend) and na (divisor) respectively 

ma=ml ; 
na=m2 ; 

% Calculate quotient and remainder 

iquot=f ix (ma/na) ; 
irem=ma-na*iquot ; 

% If remainder is not zero, reset dividend and divisor 



while irem>0 

bo3 =bol - iquot *bl ; 

bo4=bo2-iquot*b2 ; 

bol=bl ; 

bo2 =b2 ; 

bl=bo3 ; 

b2=bo4 ; 

ma=na; 

na=irem; 

iquot=fix (ma/na) ; 
irem=ma-na* iquot ; 

end 



% calculate new coefficients of ml and m2 
% redefine bol, bo2, bl and b2 

% redefine dividend and divisor 
% reapply Euclidean algorithm 
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Digital Signal Processing Solutions Products - TiMS320C54x 

JDsanHnHKBnBBBHanaBDBnmMBBnaaEaHnBHMnBHHa 

TMS320C54x DSKplus 

The 'C54x DSKplus is a low-cost design tool that gives designers a 
working knowledge of DSP code. From this foundation, designers can 
begin building complete 'C54x DSP-based systems. Priced at US $149, 
the *C54x DSKplus (part no. TMDS32000L0) is available from TI 
authorized distributors . 



The 'C54x DSKplus builds on TI's industry-leading line of low cost, 
easy-to-use DSP Starter Kit (DSK) development boards. The 
high-performance board features the TMS320C542 16-bit fixed-point 
DSP. Capable of performing 40 milli on instructions per second 
(MIPS), the 'C542 makes the 'C54x DSKplus the most powerful DSK 
development board on the market. 

Other TMS320 DSKs include the 'C2x DSK . the 'C5x DSK . and the 
floating-point 'C3x DSK . 



Key Features 

The 'C54x DSKplus includes: 

• 40 MIPS TMS320C342-based board 

• TLC320AC01 Analog Interface Circuit (AIC) 

• 'C54x DSKplus assembler, loader, Code Explorer debugger, and sample programs (3.5 M disks) 

• TMS320C54x CPU and Peripherals Reference Guide 

• TMS320C54x Algebraic Assembler Instruction Set 

• TMS320C54x Datasheet 

• TMS320C54x DSKplus User’s Guide 

• TLC320AC01 Datasheet 

• PC connector cable and universal power supply included 

• US‘5149 discount coupon toward the purchase of the 'C54x EVM 
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DSKplus Key 

Features \ 

7 


Benefits j 

i 


TMS320C542 DSP (40 
MIPS, 16-bit) 


High-performance, veiy efficient architecture requires fewer \ 

MIPS than competing DSPs to implement most algorithms 


Code Explorer debugger 
interface 


i 

/ 4 

An easy-to-use, true Windows-based interface. Supports 
symbolic debugging, breakpoints, graphical animation, variable \ 

watch windows, file I/O, algebraic/mnemonic disassembly, 
on-line help. j 

J 


Symbolic debugging (Code i 
Explorer) 


f 

Enables easy programmability by using labels for referencing 
constants, variables, matrices by name. \ 

i 


Algebraic assembler 


i 

Bypasses learning new DSP mnemonic instruction set specifics. \ 

Makes coding easier and more straight-forward. Easy one-step 
assembly and linking process. 1 

. . j 


Demo programs / 
Application code 


j 

Helps users get up-to-speed quickly j 

» 

i 


TLC320AC01 Analog 
Interface Chip 


% 

Low power dissipation, 14-bit linear resolution, programmable 
sampling rates, anti-aliasing filter, and input gain; selectable j 

auxiliary input; data read-back I 

» 

s 

j 


Socketed Programmable 
Array Logic (PAL) j 


Allows experienced designers to reprogram the PAL and change 1 
the way the host port interface works on the C54x DSKplus. 


Universal power supply & j 
cable included i 

J 


Allows for immediate use out of the box; ideal for powering 
daughter cards; filtered and regulated - thus no need for j 

on-board voltage regulation. j 



'C54x Algebraic Assembler 

The C54x DSKplus includes the algebraic assembler that speeds the initial code development process. The 
algebraic assembler does not require new users to learn a new DSP mnemonic instruction set, making 
coding easier and more direct. The assembler also utilizes a one-step assembly and linking process to 
simplify code debugging The software accomplishes this by using special directives to assemble code at an 
absolute address. 

Some extremely useful features include: 
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• In-line Assembly expression analysis allows the assembler to work when defining complex variables 
or bit locations. 

• Symbolic Debugging allows the user to reference variables by name instead of the physical address. 

• Assembling conditional blocks of assembly code using .if/.else iff.end if directives. This is especially 
helpful when you want to conditionally assemble code via a command-line argument of internal 
assembly variable. 

• Support of .sect, .bss, .usect, .text, and .data sections. 



Code Explorer Debugger 

The 'C54x DSKplus debugger was developed by GO DSP Corporation in an effort to provide the first true 
Windows-based debugger for a DSK. The Code Explorer debugger supports debugging, a new feature 
available only on the DSKplus that allows the user to specify labels for referencing constants, variables, 
and marticies by nam e Also, the debugger desktop environment is fully configurable and loaded upon 
entry into the debugger. This means that optional colors, fonts, and window sizes can be changed within 
the debugger and saved upon exiting. 

Some addi tional features of the debugger include capability of connecting files as I/O, graphical animation, 
and d ata memory viewing. The file I/O capability enables users to connect files as inputs or outputs to any 
location within your application code. Therefore you can simulate different input sequences and data 
streams without having to physically generate them. 

Graphical animation allows you to view data in a graphical format, either with time domain or frequency 
domain and in a variety of variable sizes (i.e. 8-bit signed char, 8-bit unsigned char, 16-bit, 32-bit, etc). 



Disassembly Window 






The data memory window cam be modified or 
replicated as needed. By placing the cursor inside 
the data memory window and right-clicking and 
then choosing properties, the user can change the 
title of the window, starting address and even data 
organization in the window. Valid display formats 
include 8-bit signed/unsigned char, 
signed/unsigned long, floats, and others. The page 
field can specify either Data or Program memory 
spaces. 

'C54x CPU and Peripheral 
Registers 
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OhOC 
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0x0808: 
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. _fwna->a • 


„ nofif;: 
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1 
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A 
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SI 


F~ 
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The two register windows in the 'C54x Code 
Explorer debugger are the CPU and Peripheral 
Registers. The 'C54x CPU Registers is the 
collection of registers which control the operation 
of the DSP CPU. The program counter, status 
register, and configuration registers are contained 
within this window. Notice that bit values within 
the register are brought out separately to make 
modification and monitoring easier. 

The second window is the Peripherals window. 
This window includes the registers for configuring 
the DSP peripherals like the serial ports and 
timers. Modifications to this register can be done 
by clicking on the register in the Peripheral 
Registers window'. 



| J C54X Registers \ 
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AEO = 
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B - FFFFFFFFH) 
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100b 
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AR3 * 


FFFF 


ova - 1 5a 




THH = FFFF 


AE4 = 


FFFF 


SX3 - 1 




ST0 - 1800 


AES - 


0000 


Cl 8 * 0 :i : 




ST1 = 3B00 


APS = 


FFFF 


FPCT = 0 hit 




PMST - 00 AC 


A£? - 


FFFF 


chpt * o m 




* DP = 0000 


BSC * 


FFFF 


CPI = O 




ASM =• 000C 
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AEO 
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BEAT = 0 


1ST2 * 


0 HfVi'C * 0 VX 




BHC = FFFF 
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ESA = FFFF 
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IFF - 
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i TIM = 3210 


DSP 


* oooo 


STCSE = 7FFF£t 


1 PRD - FFFF 
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- oooo 
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] TCR = 0000 
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Graphical Windows 

Graphical windows are extremely useful 
when trying to view a value of a register, 
variable, or buffer. The graphic window 
allows the user to animinate any value in 
either data or program DSP memory. 

This is accomplished by placing a 
breakpoint anywhere in the application 
code and pressing the Animation button. 

Each time the DSP reaches the 
breakpoint the graphical windows are 
updated and refreshed. 

The options window contains the 
graphics setup for the window. For 
example, the title can be changed to 
reflect the data being animated, the 
display buffer length can be changed, or 
the data read from the DSP can either be 
a single value from a list (buffer) of values in either data or program memory. Also, the sampling rate can 
be modified for correct displaying of the frequency data (EFT). The display can be viewed using 8-bit 
signed/unsigned chars, ints, long, floats, and even a log can be performed on the displayed data. 




Setting Breakpoints 
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A breakpoint can be selected by either 
double clicking on a line in the 
Disassembly window or by Selecting the 
DEBUG-BREAKPOINTS in the Pull 
Down Menu. The Pull Down Menu will 
prompt you with a menu listing all the 
available symbols in the Symbols box. 
You can either select a breakpoint from 
the list of Symbols or by entering an 
address in the Address field. 



The Breakpoint dialog box contains the 
following fields: Address, Symbols and 
Breakpoints. If the address of the desired 
breakpoint is known, simply enter the 
value in the Address field The Symbol 
field con tains the list of all the symbols in 
simply type the label name and press add 

Setting Probe Points 



Probe points allow the update of a 
particular window or the reading/wnting 
of samples from a file to occur at a 
specific point in an algorithm. This 
effectively "connects a signal probe" to 
that point in the algorithm. 

When a graph window object is created 
it assumes that it is to be updated at every 
breakpoint. However, this attribute can be 
changed and the window can be updated 
only when the program reaches the 
connected probe point. After the probe 
point is hit, and the window is updated 
execution of the program is continued 
This optimizes the display of the graph 
window and also allows you to keep a history of the signal even when the data on the DSP is not valid 

With the combination of Code Explorer's File I/O capabilities, probe points can also used to connect 
streams of data to a particular point in the DSP Code. When the probe point is reached in the algorithm, 
data is streamed from a specific memory area to file, or from the file to memory. 







y 




r 




|cSOC (fesatei* 















c ^ 




the program. If the location address of the breakpoint is labeled 



Using File I/O 




Code Explorer allows the user to stream 
data onto (or from) the target from a PC 
file. This allows the user to simulate code 
using known sample values. Note that 
this file I/O feature is not intended to 
satisfy real-time constraints. The File 
Input/Output feature uses probe points. 
When the execution of the program 
reaches a probe point, the connected 
object, whether it is a file, graph or 
memory window, is updated. Once the 
connected object is updated, execution 
continues. Using this concept, if a probe 
point is set at a specific point in the code 



and then connected to a file, file I/O functionalities can be implemented 




System Requirements 

• A 3 86, 486, or Pentium PC with a 3. 5 "disk drive 

• 4-bit parallel and/or 8-bit bidirectional parallel ports. 

• A minimum of 4Mbytes of memory 

• Color VGA monitor 

• Windows 3. 1 or Windows 95 

• ASCII editor 



How to Install 

When connecting the DSKplus to your PC, it is highly recommended you turn 
off your PC's power to make the connections below: 



1. Connect the DB25 cable (female) to the PC's Parallel port (male). 

2. Connect the DB25 cable (male) to the DSKplus board (female). 

3. Connect the power cord (NEMA cable) to the 5 volt power supply. 

4. Connect the 5-pin DIN-to-5.5mm adapter to the power supply’s 5-pin 
DIN connector. 

5. Plug the power supply power cord to the wall outlet 

6. Plug the 5.5mm connector into the power jack of the DSKplus board 

At this point the green power LED is illuminated and power is supplied to the 'C54x DSKplus board If the 
Green tpd is not illuminated, check the connections on the power supply and power cord 

Installing the software 
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The DSKplus kit includes two 3.5" floppies labeled Disk #1 and Disk #2. To 
install the software correctly, please follow the steps below: 

1. Insert Disk #1 into the 3.5" drive. 

2. From the start menn (Windows95) or the Files menu (Windows 3. 1) 
select the Run., option. Type a:\setup.exe 

3. The ins tallati on script will appear. You will be asked to select a 
destination directory. By default it will select the DSKplus directory. 

Enter the directoiy name if you would like to specify a different 
directory. 

4. When prompted, insert Disk #2 into the 3.5" floppy drive. 

5. When install ation has completed, the installation will inform you that the installa tion was 
successful. At this point a Code Explorer Group will appear. 

Starting the Debugger 



IS C:\Vr1NDOWS\SlarL-. §ii S : 
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To start the debugger, click on the icon located in the Code Explorer Group or desktop. The Code Explorer 
background and windows will appear with the Setup Box shown active. 



Select the port which is connected to the DSKplus board. If for some reason 
the port is not listed, the port address can be modified by typing in the address 
int the text box. 

As a result of selecting the correct port and proper hardware connections, the 
debu gg er will fill its windows with data and the DSKplus is now functioning. 
If for some reason the debugger responds with the error "Can't initialize 
Target DSP", follow the directions in the error box. 




T roubleshooting 



1 . Is the power on? Be sure green LED is illuminated If not, a loose power cable is ham pering your 
setup. 

2. Is the parallel part cable connection secure? In many new DSKplus boards and parallel port cables, 
substantial pressure marry be needed to connect the cables. Connect the cable to the DSKplus board 
by placing the thumb behind the DB-25 connector. Take the cable connector chassis and place 
between the index and middle fingers. Align the connectors and press the fingers together. 

3. The port selected is not being "Captured" by Windows 95. Capturing is used by Windows 95 to 
allow DOS programs access to printers. The port can be released by going into the control panel 
and selecting the printers icon. Hig hli g ht any printer and go to the File pulldown on the command 
bar. Select properties and then the Details tab. The Details tab includes a button named End 
Capture... Click on this button and select the LPT port where the DSKplus board is connected If 
the LPT is not listed, then the px»rt is not captured (select cancel) and proceed to number 4. 

4. The part selected is configured as an EPP or ECP part. The DSKplus board supports 4-bit 
unidirectional and 8-bit bidirectional parallel parts. The DSKplus does not support EPP and ECP 
parts. To check the part configuration, exit out and reboot your system. At the paint where the 
BIOS Setup routine can be selected, press the keyboard sequence to enter the BIOS (usually 
CTRL+ALT+ESC). Confirm that the parallel part is setup as '8-bit', 'bidirectional' or ’standard’ 
Specifically, not an EPP or ECP part If problems persist, run the included selftest program. 



Beyond the T C54x DSKplus 

With higher performance than any other DSK available today, the ’C54x DSKplus offers a rich 
development environment for benchmarking and evaluating code in real-time. The ’C54x DSKplus is 
designed as an easy-to-use entry into the world of high-performance fixed-point DSPs. 
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However, as your design experience grows, you may require additional functionality and expanded 
capabilities. To meet these needs, TI offers a comprehensive line of evelopment tools for the TMS320 DSPs 
that support the design process from system concept to production. 

T54x Development Tools 




© Copyright 1997 Texas Instalments Incorporated. All rights reserved 
Trademarks, Important Notice / 
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APPENDIX C 



C AND ASSEMBLY LANGUAGE CODE FOR DSP HARDWARE 



*********************************************************************** 
File: FirstApl.ASM / 

When sampling frequency is changed, need to change 

a. buffer size here 

b. A and B registers in ACOlinil.asm 

c. sampling frequency in dftsort.c 

d. buffer size in hostappl.cpp 

*********************************************************************** 
.width 80 
.length 55 

.title "FirstApp program” 

. mmregs 

.setsect n .text", 0x500,0 
.setsect "vectors”, 0x180,0 



VECTORS 



.sect "vectors" 

. copy "c : \dskplus\inits\vectors . asm" 
. text 



start : 

call AC01INIT 
pmst = #01a0h 
sp = #0ffah 
ar2 = #1200h 
*ar2 + = data(#0bh) 
imr = #280h 
intm = 0 

wait nop 

goto wait 



set up iptr 

init stack pointer. 

pointer to receive buffer at 1200h. 

store to rev buffer 

ready to rev int ' s 



XINT: 



Receive Interrupt Routine 



b = trev 
b = #0FFFCh & b 
*ar2+ = data(#0bh) ; 

tdxr = b ; 

TC = (@ar2 == #1471h) ; 

if (TC) goto restrt ; 

return enable 



load acc b with input 

store to rev buffer 
transmit the data, 
change here if fs changes 
stop if rev buffer is at 1471h 



restrt 
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ar2 = #1200h 
hpic = #Oah 



; set intm bit ...no int;'s 
; flag host task completed 



return__enable 

end isr 



. copy "c : \dskplus\f irstapp\ac01inil . asm" 

.end 

*********************************************************************** 
File: FirstApp.ASM -> First Application program for the ' C54x DSKplus 

a. buffer size here 

b. A and B registers in AC01ini2.asm 

c. sampling frequency in dftsort.c 

d. buffer size in hostapp2.cpp 

*********************************************************************** 
.width 80 
.length 55 

.title " FirstApp program" 

.mmregs 

.setsect " .text", 0x500,0 
.setsect "vectors", 0x180,0 



VECTORS 



.sect "vectors" 

.copy "c:\dskplus\inits\vectors . asm" 
. text 



start: 

call AC01INIT 
pmst = #01a0h 
sp = #0ffah 
ar2 = #1200h 
*ar2 + = data(#0bh) 
imr = #2 8 Oh 
intm = 0 

wait nop 

goto wait . 



; set up iptr 
; init stack pointer. 

; pointer to receive buffer at 1200h. 
; store to rev buffer 

; ready to rev int ' s 



Receive Interrupt Routine 



XINT: 

b = trev 
b = #0FFFCh & b 
*ar2+ = data(#0bh) 
tdxr = b 

TC = (@ar2 == #01400h) 
if (TC) goto restrt 
return enable 



load acc b with input 

store to rev buffer 
transmit the data, 
change here if fs change 
stop if rev buffer is at 1400h 
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restrt 



; set intm bit ...no int ' s 
; flag host task completed 



ar2 = #1200h 
hpic = #0ah 
return_enable 

end isr 

. copy n c : \dskplus\f irstapp\ac01ini2 . asm” 
.end 
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*********************************************************************** 



File: AC 01 INI 1 .ASM -> AC01 Initialization Routine 



*********************************************************************** 
.width 80 
.length 55 

.title " AC01 Initialization Program” 

. mmr egs 



************************************************************************ 



* Certain AC01 registers can be initialized using a conditional assembly 

* constant. By setting the constant REGISTER to the appropriate value, 

* the assembler will either include initialization for certain registers 

* or ignore register initialization. 



* The constant REGISTER should be set to include the following AC01 

* register: 

* 

* REGISTER (binary) = 



* 

* 

* 

* 

★ 

* 

★ 

* 

* 

★ 

* 

★ 



0000 


0000 


0000 


0001 


- > 


initialize 


Register 


1 


(A Register) 


0000 


0000 


0000 


0010 


- > 


initialize 


Register 


2 


(B Register) 


0000 


0000 


0000 


0100 


- > 


initialize 


Register 


3 


(A' Register) 


0000 


0000 


0000 


1000 


- > 


initialize 


Register 


4 


(Amplifier Gain- 
Select) 


0000 


0000 


0001 


0000 


- > 


initialize 


Register 


5 


(Analog 

Conf igurat ion) 


0000 


0000 


0010 


0000 


- > 


initialize 


Register 


6 


(Digital 

Configuration) 


0000 


0000 


0100 


0000 


- > 


initialize 


Register 


7 


(Frame -Sync Delay) 


0000 


0000 


1000 


0000 


- > 


initialize 


Register 


8 


(Fr am -Sync number) 



* Any combination of registers can be initialized by adding the binary 

* number to the REGISTER constant. For example to initalize Registers 4 

* and 5 , REGISTER = I8h. Upon assembly, only code for register 4 & 5 

* initialization is included in the AC01INIT module. When called the 

* module will load REG4 and REGS values into internal AC01 registers. 



* 



* 

* Register 4 is always loaded to get a 6db input gain. This sets full- 

* scale to 3v (p-p input) due to the single-ended AC01 configuration. 

* 

* 



REGISTER 


.set 


Obh 


; Powerup 


default values 


REG1 


.set 


lfeh 


. * 
t 


H2h 


REG2 


. set 


2lfh 


. * 
t 


212h 


REG3 


. set 


3 0 Oh 


t 


300h 


REG4 


. set 


4 0dh 


. * 


405h 


REG5 


. set 


501h 


/ 


501h 


REGS 


. set 


SOOh 


t 


SOOh 
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REG7 

REG8 



.set 700h 
.set 801h 



700h 

801h 



AC01INIT: 

xf = 0 
intm = 1 
ter = #10h 
imr = #2 8 Oh 
tspc = #0008h 
tdxr = #0h 
tspc = #00c8h 
xf = 1 

; Register 



reset acOl 

disable all int service routines 
stop timer 

wakeup from idle when TDM Xmt int 
stop TDM serial port 
send 0 as first xmit word 
reset and start TDM serial port 
release acOl from reset 

init's 



.eval REGISTER & lh, SELECT 

.if SELECT = lh 

a = #REG1 

call REQ2 

.endif 



if REG1 then include this source 

load Acc A with REG1 value 
Call REQ2 subroutine 



.eval REGISTER & 2h, SELECT ; if REG2 then include this source 
.if SELECT = 2h 
a = #REG2 
call REQ2 

.endif 



.eval REGISTER & 4h, SELECT ; if REG3 then include this source 

.if SELECT = 4h 

a = #REG3 

call REQ2 

.endif 

.eval REGISTER & 8h, SELECT ; if REG4 then include this source' 

.if SELECT = 8h 

a = #REG4 

call REQ2 

.endif 



.eval REGISTER & 10h, SELECT ; if REG5 then include this source 

.if SELECT = lOh 

a = #REG5 

call REQ2 

.endif 

.eval REGISTER & 2 Oh, SELECT ; if REG6 then include this source 

.if SELECT = 20h 

a = #REG6 

call REQ2 

.endif 
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. eval REGISTER & 4 Oh, SELECT ; if REG7 then include this source 

. if SELECT = 40h 

a = #REG7 

call REQ2 

.endif 

.eval REGISTER & 80h, SELECT ; if REG8 then include this source 

.if SELECT = 8 Oh 

a = #REG8 

call REQ2 

. endif 

return 



REQ2 



ifr = #080h 
tdxr = #03h 



; clear flag from IFR 
; request secondary when AC01 starts 



idle(l) 
tdxr = a 
ifr = #080h 



; wait for primary to xmit 
; send register value to serial port 
; clear flag from IFR 



idle (1) 
tdxr = #0h 
ifr = #080h 
idle (1) 



wait for secondary to xmit 

send neutral state in case last init 

clear flag from IFR 

wait for neutral state to xmit 

return from subroutine 



return 
. end 
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.*********************************************************************** 



; File: AC01INI2 .ASM -> AC01 Initialization Routine 



.**************************************************************★******** 



.width 80 
.length 55 

.title "AC01 Initialization Program" 
.mmregs 






* Certain AC01 registers can be initialized using a conditional assembly 

* constant. By setting the constant REGISTER to the appropriate value, 

* the assembler will either include initialization for certain registers 

* or ignore register initialization. 



* The constant REGISTER should be set to include the following AC01 

* register: 

* 

* REGISTER (binary) = 

* 



* 

* 

★ 

* 

* 

* 

* 

★ 

★ 

★ 

* 

* 



0000 


0000 


0000 


0001 


- > 


initialize 


Register 


1 


(A Register) 


0000 


0000 


0000 


0010 


- > 


initialize 


Register 


2 


(B Register) 


0000 


0000 


0000 


0100 


- > 


initialize 


Register 


3 


(A’ Register) 


0000 


0000 


0000 


1000 


- > 


initialize 


Register 


4 


(Amplifier Gain- 
Select) 


0000 


0000 


0001 


0000 


-> 


initialize 


Register 


5 


(Analog 

Configuration) 


0000 


0000 


0010 


0000 


- > 


initialize 


Register 


6 


(Digital 

Configuration) 


0000 


0000 


0100 


0000 


- > 


initialize 


Register 


7 


(Frame-Sync Delay) 


0000 


0000 


1000 


0000 


-> 


initialize 


Register 


8 


(Fram-Sync number) 



* Any combination of registers can be initialized by adding the binary 

* number to the REGISTER constant. For example to initalize Registers 4 

* and 5, REGISTER = 18h. Upon assembly, only code for register 4 & 5 

* initialization is included in the AC01INIT module. When called the 

* module will load REG4 and REG5 values into internal AC01 registers. 



* 



* Register 4 is always loaded to get a 6db input gain. This sets full- 

* scale to 3v (p-p input) due to the single-ended AC01 configuration. 



* 



REGISTER 


.set 


Obh 


REG1 


. set 


If eh 


REG2 


. set 


23ch 


REG3 


.set 


300h 


REG4 


.set 


4 0dh 


REGS 


. set 


501h 


REG6 


. set 


600h 



Powerup default values: 
112h 
2 12h 
300h 
405h 
501h 
600h 
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REG7 .set 700h 

REG8 .set 801h 



700h 
8 Olh 



ACO UNIT : 

xf = 0 
intm = 1 
ter = #10h 
imr = #2 8 Oh 
tspc = #0008h 
tdxr = #0h 
tspc = #00c8h 
xf = 1 

; Register 



reset acOl 

disable all int service routines 
stop timer 

wakeup from idle when TDM Xmt int 
stop TDM serial port 
send 0 as first xmit word 
reset and start TDM serial port 
release acOl from reset 

nit’s 



.eval REGISTER & lh, SELECT ; if REG1 then include this source 
.if SELECT = lh 

a = #REG1 ; load Acc A with REG1 value 

call REQ2 ; Call REQ2 subroutine 

.endif 

.eval REGISTER & 2h, SELECT ; if REG2 then include this source 
.if SELECT - 2h 
a = #REG2 
call REQ2 



.endif 



.eval REGISTER & 4h, SELECT ; if REG3 then include this source 

.if SELECT = 4h 

a = #REG3 

call REQ2 

.endif 

.eval REGISTER & 8h, SELECT ; if REG4 then include this source 

.if SELECT = 8h 

a = #REG4 

call REQ2 

.endif 

.eval REGISTER & lOh, SELECT ; if REG5 then include this source 

.if SELECT = lOh 

a = #REG5 

call REQ2 

.endif 



.eval REGISTER & 20h, SELECT ; if REG6 then include this source 
.if SELECT = 20h 
a = #REG6 
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call REQ2 
. endif 

. eval REGISTER & 40h, SELECT ; if REG7 then include this source 

.if SELECT = 4 Oh 

a = #REG7 

call REQ2 

.endif 

.eval REGISTER & 80h, SELECT ; if REG8 then include this source 

.if SELECT = 8 Oh 

a = #REG8 

call REQ2 

. endif 

return 



REQ2 



if r = #080h 
tdxr = #03h 



; clear flag from IFR 
; request secondary when AC01 starts 



idle(l) 
tdxr = a 
ifr = #080h 



; wait for primary to xmit 
; send register value to serial port 
; clear flag from IFR 



idle(l) 
tdxr = #0h 
ifr = #080h 
idle (1) 



wait for secondary to xmit 

send neutral state in case last init 

clear flag from IFR 

wait for neutral state to xmit 

return from subroutine 



return 

.end 
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y**********************************************************************^/ 

/* 

/* File: H0STAPP1 . CPP Source code for host application 
/* 

# include <HI54X.H> 

#include <stdio . h> 

#include <conio .h> 

#include <stdlib .h> 

extern int datareg [] , statreg [] , ctrlreg [] ; 
extern int pport , portmode , Readdelay; 



void main (void) 

{ 



FILE *fp ; 

if ( (fp=f open ("datal.dat" , "w") ) = 



=NULL) 



{ 



clrscr ( ) ; 

pr int f ( "Cannot open file 
exit (0) ; 



An") 



/* Open file */ 



portmode=0; 

Readdelay = 20; 
clrscr () ; 

if ( (pport =locat export () ) >= 5) { 

printfC'No connection\n 
backout () ; 
exit (0) ; } 

else{ } 

_setcursortype (_NOCURSOR) ; 
set_latch (1,1) ; 
int word =0, col=0; 



/* 

/* 

/* 

/* 



Col=0 ; 

gotoxy (1, 1) ; 

send_word (0x0808 , C_SEND) 
HINT (10000) ; 

send word (0x1200, A SEND) 



4 -bit mode */ 

In case host slow*/ 
Clear the screen */ 
Find the port. */ 
/* If no connection 
/* then leave board 
/* in known state 



/* Hide text cursor * 
/* Keep DSP running * 
/* and bring PAL out* 
/* out of Tri-state * 



for (int buf=0 ; buf < 



{ 



/* go to home * 

/* Clear the HINT * 
/* Wait for nxt HINT* 
/* Goto 0x46 entries* 
/* before buffer * 
0x271; buf++) 

/* change here if fs change* 



word = read_word(D_READ); 
print f ( ”%4 .4x ”, word) ; 
fprintf(fp, "%d\n", word) ; 
if (col >= 13) { 

COl=0; 

print f ( "\n fl ) ; } 
else{col++; } 

} 

_setcursortype ( _NORMAL CURS OR ) ; 



/* 

/* 

/* 

/* 



Read word from pp* 
Print it to scr * 
Output to file * 
in 14 columns * 



/* Ret normal cursor* 



/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



/ 
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} 



fclose (fp) ; 
backout ( ) ; 

exit (0) ; 



/* Close file */ 

/* Leave board in */ 
/* known state */ 
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^********************************************************************** j 
/* 

/* File: H0STAPP2.CPP Source code for host application 
/* 

^********************************************************************** j 



# include <HI54X.H> 
#include <stdio.h> 
#include <conio . h> 
#include <stdlib.h> 



extern int datareg[] , statreg[] , ctrlregU ; 
extern int pport, portmode, Readdelay ; 



void main (void) 

{ 



FILE *fp ; 

if ( (fp=f open ("data2.dat", "w") )==NULL) 



/* Open file */ 



clrscr ( ) ; 

printf ( "Cannot open file .\n") 
exit (0) ; 



portmode=0 ; 

Readdelay = 20; 
clrscr ( ) ; 

if ( (pport=locate_port ( ) ) 



/* 
/* 
/* 
5) {/* 



4 -bit mode */ 

In case host slow*/ 
Clear the screen */ 
Find the port. */ 



printf ("No connection\n" ) ; /* 



backout ( ) ; 
exit (0) ; } 

else{ } 

_set cursor type (_N0CURS0R) ; 
set_latch(l, 1) ; 
int word =0, col=0; 



/* 

/* 

/* 

/* 

/* 

/* 



If no connection */ 
then leave board */ 
in known state */ 

Hide text cursor */ 
Keep DSP running */ 
and bring PAL out*/ 
out of Tri-state */ 



col=0 ; 

gotoxy (1, 1) ; 
send_word (0x0808 , 
HINT (10000) ; 
send word (0x1200, 



C_SEND) 
A SEND) 



/* 

/* 

/* 

/* 

/* 



for (int buf=0 



{ 



go to home */ 

Clear the HINT */ 
Wait for nxt HINT*/ 
Goto 0x46 entries*/ 
before buffer */ 
buf < 0x200; buf++) 

/* Change here if fs is changed*/ 



word = read_word (D_READ) ; 
printf ( " %4 . 4x " , word) ; 
fprintf(fp, "%d\n" , word) 
if (col >= 13) { 

COl=0 ; 

printf ( "\n" ) ; } 
else{col++; } 

} 



/* 

/* 

/* 

/* 



Read word from pp*/ 
Print it to scr */ 
Output to file */ 
in 14 columns */ 
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} 



_setcursortype (_NORMALCURSOR) ; 
fclose (fp) ; 
backout { ) ; 
exit (0) ; 



/* Ret normal cursor*/ 
/* Close file */ 

/* Leave board in * / 
/* known state */ 
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^*ic*icic*ic1c*icic1c***icic-k-k-kic-kir-k1c-k1c-k-k1c*ic-k-kic***ic-k*-kic-k1c-k-k-k*1c-k-k*ic-kic*-k-k*ic-k*-k-k’k*-k-k* j 

/* 

/* File: MAIN.C Source code for main program 
/* 

j*-kic*ic**ic********icicic1c****ic***icic1c-k1cicicic1c*-k*ic-kicic*ic-kic**-k***1c1c*icic1c*-kic1c**1c1cic* j 



# inc lude "math.h" 
#include " s tddef . h" 
#include "stdlib.h" 
# inc lude "stdio . h" 

# inc lude " conio . h " 



#define fsl 625 
#define fs2 324 

main ( ) 

{ 

FILE *fpl, * f p2 ; 

int dl, d2, xl [fsl] , x2 [fs2] ; 

int k, out, freq; 

int glde(int k) ; 

double xrolffsl], xiol [fsl] , xol [fsl] ; 

float pi = 3.1415926, tpi; 

int n, u; 

int i, dftl, nl; 

double maxi , max2 ; 

double xro2 [fs2] , xio2 [fs2] , xo2 [fs2] ; 
int j, dft2, n2; 

if ( (fpl=fopen( n data5.dat” , "rt") )==NULL) /* Open file */ 

{ 

clrscr () ; 

pr int f ( "Cannot open file .\n ,T ); 
exit (0) ; 

} 

if ( (fp2=fopen ( f, data6 .dat " , "rt " ) ) ==NULL) /* Open file */ 

{ 

clrscr () ; 

printf ("Cannot open file . \n") ; 
exit (0) ; 

} 

for (n=0 ;n<f si ;n++) 

{ 

fscanf(fpl, "%e ", &dl) ; 
xl [n] =dl ; 

} 

for (n=0 ;n<f s2 ;n++) 

{ 

fscanf(fp2, "%e ", &d2) ; 
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x2 [n] =d2 ; 

} 

fclose (fpl) ; 
f close (fp2) ; 

tpi=2*pi ; 

for (u=0 ;u<f si ;U++) 

{ 

xrol [u] = 0.0; 
xiol [u] = 0.0; 
for (n=0 ;n<fsl ;n++) 

{ 

/*-- Xr [u] = (1/fsl) sum (xr [n] . cos (2PI .u .n/f si) } --*/ 
xrol [u] = xrol [u] + xl [n] *cos ( tpi*u*n/f si ) ; 
/*-- Xi [u] = - (1/fsl) sum xr [n] . sin (2PI .u .n/f si) --*/ 

xiol [u] = xiol [u] - xl [n] *sin (tpi*u*n/f si) ; 

} 

xrol [u] =xrol [u] / f si ; 
xiol [u] =xiol [u] / f si ; 

xol [u] =sqrt (xrol [u] *xrol [u] +xiol [u] *xiol [u] ) ; 

} 

dftl=0 ; 
nl=f sl/2+1 ; 
maxl=xol [0] ; 
for (i=l ; i<nl ; i++) 

{ 

if (xol [i] > maxi) 

{ 

df tl = i ; 
maxl=xol [i] ; 

} 

} 

for (u=0 ; u<f s2 ;U++) 

{ 

xro2 [u] =0.0; 
xio2 [u] =0.0; 
for (n=0 ; n<f s2 ; n++) 

{ 

xro2 [u] = xro2 [u] + x2 [n] *cos ( tpi*u*n/f s2) ; 
xio2 [u] = xio2 [u] - x2 [n] *sin (tpi*u*n/f s2) ; 

} 

xro2 [u] =xro2 [u] /fs2; 
xio2 [u] =xio2 [u] /f s2 ; 

xo2 [u] =sqrt (xro2 [u] *xro2 [u] +xio2 [u] *xio2 [u] ) ; 



df t2=0 ; 
n2=f s2/2 ; 
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max2=xo2 [ 0 ] ; 
for ( j =1 ; j <n2 ; j++) 

{ 

if (xo2 [ j ] > max2) 

{ 

df t2=j ; 
max2=xo2 [ j ] ; 

} 

} 

k=df t2-df tl ; 

out=glde (k) ; 

freq=abs (dftl+out*fsl) ; 

printf( T ’The frequency is %d", freq) ; 

for ( ; ; ) ; 

} 

glde (k) 

{ 

float mult / bl / b2 , mtest / mdl / md2 / mx; 

/* This section solves the linear diophantine equation fsl*bl + fs2*b2 
na where fsl and fs2 are the sampling frequencies 

and na is the greatest common divisor and returns the value bl, b2 and 
na. fsl and fs2 are assumed positive */ 

float bol , bo2 ,ma, na , irem,bo3 ,bo4 ; 
int iquot ; 



bol=l ; 
bo2=0 ; 
bl = 0 ; 
b2=l ; 

/* Place fsl and fs2 in ma (dividend) and na (divisor) respectively */ 

ma=fsl; 
na=fs2 ; 

/* Calculate quotient and remainder */ 

iquot=ma/na; 
irem=ma-na* iquot ; 



/* If remainder is not zero, reset dividend and divisor */ 

while (irem>0) 

{ 

bo3=bol- iquot *bl ; 
bo4=bo2- iquot *b2 ; 
bol-bl; 
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/* calculate new coefficients */ 

/* redefine bol, bo2, bl and b2 */ 



} 



bo2=b2 ; 

bl=bo3 ; 

b2=bo4 ; 

ma=na; 

na=irem; 

iquot=ma/na; 

irem=ma-na*iquot ; 



/* redefine dividend and divisor */ 
/* reapply Euclidean algorithm */ 



/* To check whether the equation is solvable, na must be a factor of k 
for the equation to be solvable. */ 



mult=k/na; 
if ( (k-mult*na) ==0) 
{ 

bl=bl*mult ; 

b2=b2*mult ; 
mtest=bl; 

mdl=fsl/na; 
md2=fs2/na; 
mx=bl; 
mx=mx+md2 ; 



/* Equation is solvable */ 



/* These new values solve the */ 

/* diophantine equation */ 

/* To check whether bl and b2 */ 

/* are the least values that */ 



/* satisfies the diophantine equation */ 



while ( (abs (mx) -abs (bl) ) < 0) 

{ 

bl=mx; 
b2 =b2 -mdl ; 
mx=mx+md2 ; 

} 



if 

{ 



} 

} 

return {(int)bl); 

} 



( (mtest-bl) — = 0 ) 

mx=bl ; 
mx=mx-md2 ; 

while ( (abs (mx) -abs (bl) ) <0) 

{ 

bl=mx; 
b2 =b2+mdl ; 
mx=mx-md2 ; 

} 
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